An efficient and long-time accurate third-order algorithm for the Stokes-Darcy system

A third-order in time numerical IMEX-type algorithm for the StokesDarcy system for flows in fluid saturated karst aquifers is proposed and analyzed. A novel third-order Adams-Moulton scheme is used for the discretization of the dissipative term whereas a third-order explicit Adams-Bashforth scheme is used for the time discretization of the interface term that couples the Stokes and Darcy components. The scheme is efficient in the sense that one needs to solve, at each time step, decoupled Stokes and Darcy problems. Therefore, legacy Stokes and Darcy solvers can be applied in parallel.The scheme is also unconditionally stable and, with a mild time-step restriction, long-time accurate in the sense that the error is bounded uniformly in time. Numerical experiments are used to illustrate the theoretical results. To the authors’ knowledge, the novel algorithm is the first third-order accurate numerical scheme for the Stokes-Darcy system possessing its favorable efficiency, stability, and accuracy properties. This work is support in part by a grant from the NSF. Wenbin Chen School of Mathematical Sciences, Fudan University E-mail: [email protected] Max Gunzburger Department of Scientific Computing, Florida State University E-mail: [email protected] Dong Sun Department of Mathematics, Florida State University E-mail: [email protected] Xiaoming Wang Department of Mathematics, Florida State University E-mail: [email protected] Corresponding author 2 Wenbin Chen et al.